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1. Introduction 

The high energy scattering in the perturbative QCD has been extensively studied during the last two decades. Special 
interest is drawn to the situation where the scattering amplitude unitarizes, which is called a saturation regime ^ ^ . 
In this regime the density of partons becomes so high that one cannot ignore their mutual interactions anymore. 
A very convenient framework for studying parton saturation is the so-called dipole approach where the size 
of the elementary degree of freedom, namely, dipole is assumed to be much smaller than 1/Aqcd justifying the 
use of perturbative QCD. A lot of observables were calculated using this extremely convenient formalism. In the 
present paper we consider the Abramovsky-Gribov-Kancheli (AGK) jsj cutting rules in the framework of the dipole 
approach. The AGK cutting rules allow us to expand the dipole approach for the calculation of the scattering 
elastic amplitude (total cross section) to the consideration of exclusive processes such as diffractive production 
and different correlations in multiparticle production processes. They lead also to a better understanding of the 
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fct-factorization in the region of low Xsjorken where the new momentum scale: saturation momentum, makes this 
factorization questionable. Being such an important tool, the AGK cutting rules have been discussed in QCD (see 
Refs. 0), but the recent study of the inclusive cross section in QCD has drawn a new attention to these rules. The 
project was started by Kovchegov [|j and Kovchegov and Tuchin ||] for the single inclusive cross section and it was 
extended to the case of the double gluon production by Kovchegov and Jalilian-Marian . The main result of this 
study could be summarized as follows: for the single inclusive cross sections the /ct-factorization has been proved 
which leads to a plausible validity of the AGK cutting rules while in the case of the double inclusive cross section 
the explicit violation of the AGK cutting rules has been found. 

The result of Refs. [| | for the single inclusive cross sections was confirmed afterwards by Braun [ pT| using the 
reggeized gluon technique, as well as, by Marquet [|l2| and Kovner and Lublinsky [ p^ in the Wilson lines formalism. 

The result of Ref. was strongly questioned by Braun who claimed that one of the gluons (upper) is 
necessarily emitted from the vertex and such a contribution cannot represent a genuine violation of the AGK cutting 
rules, because the original derivation of the AGK rules was based on the assumption that there are no emissions from 
vertices. 

Motivated by this discrepancy in the results, we revisited first the single inclusive case and found that the result 
of Kovchegov and Tuchin is correct. In Section]^ we discuss the key problem of real- virtual cancellation which is the 
essential question for approaching the AGK cutting rules problem in QCD. Indeed, the key question for everybody 
and the nightmare of everybody who is interested in the AGK cutting rules is whether the set of diagrams that 
contribute to the total cross section, and to the inelastic production is the same. In this section as well as in the 
sections we complete the analysis, started by Chen and Mueller and by Kovchegov and Tuchin[^ of the 
set of diagrams that are responsible for the total and inelastic cross sections. 

In Section js] and Section ^ we repeat the Kovchegov and Tuchin analysis of the single inclusive production but 
introducing a new function M {jk\ik) in which it is easy to separate elastic and inelastic interaction. We obtain a new 
non-linear evolution equation and solve it. Function M {jk\ik) describes the non-diagonal process of scattering of 
dipole with size rj^ = \xj — Xk \ off the target with transition of this dipole to the dipole of the size r^f^ = \xi — Xh\. At 
first sight such an amplitude contradicts the key idea of the dipole approach: dipoles are correct degrees of freedom at 
high energy or, in other words, the interaction matrix is diagonal with respect to the sizes of the interacting dipoles. 
We show by solving the evolution equation for M {jk\ik), that the dipole approach survives this test of existence the 
non-diagonal amplitudes at high energy. 

Starting from Section |^ wc consider the AGK cutting rules. The AGK cutting rules give us the relation between 
the total cross section at high energy and the processes of multiparticle production. The main idea stems from the 
unitarity constraint for the BFKL Pomeron fill] which describes the high energy scattering amplitude in the leading 
log (l/xBjorken) approximation of perturbative QCD. The unitarity reads as 

2N{Y;x,y) - \N iY;x,y) \^ + iY;x,y) (1.1) 

where Y = (l/x) and (x, y) are the coordinates of the incoming dipole; N (Y; x, y) is the imaginary part of the elastic 
amplitude and the first term describes the elastic scattering (assuming that the real part of the amplitude is small at 
high energy) while the second stands for the contribution of all inelastic processes. In the leading log {1/ XBjorken) 



approximation the elastic contribution can be neglected and for the BFKL Pomeron Eq. (1.1) can be reduced to the 
form (see Fig. ^: 

2N^^^^{Y;x,y) = Gf,r'{Y;x,y) (1.2) 
In what follows we call Gin iX'^ v) ^ cut Pomeron while j^bfkl j-y . ^-j .^jj^j called a Pomeron or uncut Pomeron. 
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Figure 1: 

ladder. 



Cut Pomeron 

The definition of cut Pomeron tiirougii tiie BFKL 



In this paper we prove the AGK cutting rules for 
the case of the totally inclusive processes (see Fig. ||). 
We show that the same set of the diagrams determine 
the total cross section and the processes of the inelastic 
production. Therefore, each inelastic contribution for 
the exchange of two BFKL Pomerons can be calculated 
in terms of the elastic amplitude with the coefficients 
shown in Fig. ^ The situation changes drastically for 
the case when we measure one extra particle (see Fig. |[ 
In this case the inelastic production stems from different 
set of the diagrams than the total cross section and 
we do not have the simple expression for each inelastic 



processes through the elastic amplitude as well as we do not find the simple relations between different processes. 
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Figure 2: The AGK cutting rules for the total inclusive processes. Cut Pomeron is defined in Fig. hi and in Eq. (1.2) 





Figure 3: Violation of AGK cutting rules for the processes with additional emission of one particle (gluon). Cut Pomeron is 
defined in Fig. |l| and in Eq. (1.2). 



In the Appendices we adduce some calculations relevant to our discussion. 



2. Real- virtual cancellations 



In this section we briefly review the main result of Ref. with a special emphasis on so-called real- virtual cancel- 
lations. These cancellations play an important role in the derivation of both BFKL flB] and BK |17, 1^ equations in 



the colour dipole model. Let us consider a colourless onium state with one extra soft gluon emission. The soft gluon 
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Figure 4: All possible emissions of a soft gluon in the onium state, diagrams B* and C* are not shown. 



can be emitted either from a quark or from a antiquark lines. For simplicity we consider here only the emissions from 
the antiquark line. We assign the transverse coordinates xi, xq and X2 to the quark, antiquark and soft gluon lines, 
respectively. The system interacts with the target by instantaneous interaction via Coulomb gluon exchange. We 
use an eikonal approximation allowing for multiple interaction with the target. It was shown by Kovchegov Q that 
there are no soft gluon emissions during the interaction time so that the eikonal rescattering can be regarded as an 
instantaneous as well. To avoid any question regarding the use of the eikonal approximation as the initial condition 
for the evolution at low energy, we consider the interaction with heavy nuclei for which the eikonal formula, as well 
as, the non- linear evolution equation in the mean field approximation (Balitsky-Kovchegov (BK) equation [0 
are proved. We work in Light Cone Perturbation Theory (LCPT)jl9| with hght-cone gauge and denote by r = the 
interaction time while the detector that measures the particles in the final state is placed at r = oo. 

For simplicity, we consider only soft gluon emissions from the antiquark with coordinate xq. The relevant 
diagrams are shown in Fig. ^ 

In diagram A of Fig. ^ the soft gluon at X2 is emitted before the interaction time in both amplitude and conjugate 
amplitude. In diagrams R and R* the soft gluon is emitted and absorbed before the interaction time, these diagrams 
give reggeization term of the BFKL and BK equations in the dipole model. Only diagrams A, R and R* are present 
in the BFKL|]l^ and BK ^ equations, all other diagrams are canceled as we will show shortly. This cancellation 
is called real-virtual cancellation, where any real emission (absorption) of the soft gluon after time r = is canceled 
by its virtual counterpart (s). Namely, let us consider diagrams B and C where the soft gluon is emitted before the 
interaction time in the amplitude, but not in the conjugate amplitude. The emitted gluon be either present or not 
in the final state at r = oo. As it was shown in Ref. ( see also Appendix A for more details) those diagrams 
differ only by a minus sign and thus cancel each other in the total cross section (same for B* and C*, not shown 
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here). Similar cancellation happens also to D, E and but one should note that E and F have a factor of ^ w.r.t. 
diagram D due to a light-cone time ordered integral as shown in Appendix of Ref. . 

Thus we are left with diagrams R and R* . To this we should add also diagrams with the soft giuon emission 

from the quark line at xi, which translates into the BFKL kernel in the dipole model, namely, f ^ — ^ j = f^"^ . 

In Section ^ we present the derivation of BK equation in more general case and show how this reduces to 
well-known equation 



(97V(10) as [ ,2 X 



df'x2 2 2 {^(12) + A^(02) - A^(IO) - jV(12)jV(02)} (2.1) 

with A^(IO) = linAei{xi,xo;Y) being the imaginary part ^ of the elastic scattering amplitude of a colourless dipole 
with quark and antiquark coordinates at xi and xq (xiq = xi ~ xq), respectively. 



3. Inclusive one gluon production: no evolution included 

In this section we follow the lines of Ref. in the derivation of the single gluon inclusive production cross section in 
DIS with no evolution included. However, in our derivation we treat separately all the contributions that sum into 
the final result of Ref. The reason for that is our goal to include properly the evolution in the inclusive cross 
section for one gluon production as will be clarified in Section ^. 

The general expression for the gluon production cross section in DIS is given by 

= A / d?x,, dz\r -'^^{xouz) r r (3.1) 



<Pk dy 2Tr^ J ' d^k dy 

where d<7'^^^{xi, xo)/d'^k dy is the gluon production cross section for the scattering of a colour dipole with quark 
coordinate xi and antiquark coordinate xq on the target, and ip'^ ^'^''{xqi, z) is the well-known [ pO| wave function of 
the splitting of the virtual photon in DIS into qq pair with of a transverse size xiq = xi — Xq and a fraction z of the 
longitudinal momentum of the photon carried by the quark. 

We want to point out that we limit our discussion to a case where the produced gluon is the hardest gluon 
emitted in qq system. 

We start with selecting only those diagrams from Fig. ^, which have a real gluon emission ( the soft gluon appears 
at the cut at r = 00). It is easy to see that we are left with A, B (B*) and D. In the inclusive cross section the 
transverse momentum k of the soft gluon is kept fixed and thus transverse coordinates of the produced gluon are 
different in the amplitude {X2) and the conjugate amplitude {x2')- Before writing the expression for the inclusive cross 
section with one gluon production we define an object which will play a central role in our further derivations. The 
function M(12|34) is defined as an unintegrated over the impact parameter cross section of the scattering of a dipole 
with coordinates of quark (antiquark) being xi (X2) in the amplitude and X3 (x^) in the conjugate amplitude. One 
should keep in mind that the cross section M (12|34) is a function of rapidity though it is not reflected in our notation. 
The explicit expression for Af(12|34) = Afo(12|34) in the case of interaction with a nucleus without evolution is given 
in Appendix In fact, we have three of them, for each of the diagrams A, B {B*) and D in Fig. ^. In Section |^ 
we show that all three functions are described by the same evolution equation and obey a generalized form of the 
optical theorem for the case of the dipole having different sizes in the amplitude and the conjugate amplitude. 



^For simplicity, we omit the argument Y = log{l/xgjarken) where xgjarkan is the fraction of energy carried by the dipole. 
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We are ready to write the expression for the inclusive cross section in terms of M"*, and as follows ^ 



do-99^^<^^(xi,xo) _ as 1 r 
^Pkdy ~ 2^ (27r)2 J ' 



\ (^12 


X02\ 


[ ^12' 


_ ^02' \ 
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^02 / 


[-I2' 


^02' / 



+Mo(12|12')Mo(20|2'0) - Afo(12|12') {Afo(20) + iVo(2'0)} - A/o(20|2'0) {Afo(12) + Afo{12')} + A'o(20)A^o(12') + JVo(12)A'o(20')) (3.3) 

'xi2 _ xo2\ L _ |Afo(i2|lo) (1 - Afo(20)) + 7Vo(20)7Vo(10)} 

^^12 ^02 / V ^'02' / 



(3.4) 
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(3.5) 



3^12' a:o2' 



^^12' 



''02' 



- ) {Mo(10|2'0) (1 - Afo(12')) + Afo(10)Afo(12')} 



12' 



''12 ■ 
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''02' 



"-02 



Afo(10|10)-ipip^ + pi^l7V|(10)} 

\^Xj2 J-02' 02 ■''12' / 

where Nq is the initial condition of BK equation and M{){ik\il) is given in Appendix ^ 



(3.6) 



In Eq. ( 3.2 ) we account for a separate rescattering of dipoles " 12" (" 12"' ) and " 20" (" 2'0" ) coming from diagram A. 

In Eq. (3.3) we include the contributions where both of the dipoles are rescattered at the same time in diagram 
A. The first term describes the case where the two dipoles are rescattered both elastically and inelastically. The 
second and the third terms represent both elastic and inelastic rescattering of one dipolc and only elastic (cither 
amplitude or conjugate amplitude) rescattering of the other dipole. The fourth and fifth terms describe only elastic 
rescattering of the two dipoles. One should note that terms Nq{2Q)Nq{\2') + A'^o(12)iVo(20') are not included into 
Mo(12|12')Mo(20|2'0) since the elastic part of the cross section Mo(12|12') by the definition starts with the lowest 
order of A'^o(12)iVo(12') (there are no disconnected diagrams), and thus the terms iVo(20)7Vo(12') + A'^o(12)A'^o(20') 
should be added separately. 

In Eq. ( ^.4[ ) we have contributions from diagram B{B*) which are not symmetrical in the Kernel. Namely, some 

diagrams are prohibited by the definition of the cross section Mo(12|10). As a simple example, we try to include 

contributions which do not enter Eq. (3.4). We consider ^^^^^Mq{12\1Q) as shown in Fig. pi It is immediately clear 

^ — ' ^12 ^12' 

that Mo(12|10) cannot multiply such a dipole splitting, since a dipole "12" ("10") is not present at the interaction 
time T = in the conjugate amplitude with such a dipole splitting ^^iiai^ other words, the dipole formed by the 
lower quark loop scatters only elastically. 



The first term in Eq. (3.4) includes elastic and inelastic rescattering of only one dipole, and also elastic rescattering 
of the other dipole. The term 7Vo(20)A^o(10) corresponds to the elastic scattering of the dipole "20" in the amplitude 
and the elastic scattering of the dipole " 10" in the conjugate amplitude, as it was mentioned before this term is not 
present in M(20|10)7V(20) by the definition of M(20|10) due to the absence of the disconnected diagrams. Same for 



all other terms in Eq. (3.4) and Eq. (3.5) 



In Eq. (3.2) and everywhere in this paper we use notation kx instead of ■ x^. We hope that this will not cause difficulty in 



understanding. 
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T= T=C0 T=0 



Figure 5: An example of a diagram that leads to asymmetric Kernels in Eq. ( p.4| ) and Eq. (3.5) 



Finally, in Eq. (3.6) we account for the contribution from diagram D, where we have both elastic and inelastic 
scatterings of the initial dipole " 10" . This can be easily understood from all D-type diagrams, where the measured 
gluon is emitted from quark and/or antiquark after the interaction time r = 0. In the case of "crossed" dipole 
splittings ^2.^021 etc. the inelastic rescattering is suppressed in large Nc, that is the reason why the last term in 

Eq. ( |3.6|) has iVg (10) instead of A/(10|10)(which also includes the inelastic part). 

As the last step in our derivation we have to substitute the expressions for the total cross section Mq and the elastic 
amplitude A^o by their eikonal formulae which are found in Appendix |^. After some tedious, but straightforward 
algebra we obtain 



(Pkdy ~ 2^ (27r)2 



:2a;2d2^2,e-*^-(-^--^') (l + e-2-^2'Q?/4 _ ^-^-loQl/^ _ e-2-^'oQ^/4) (3.7) 



^12 ^02' 



(3.8) 



X02 



(,Q2 J.]^2 



e ^■^02''- 



(3.9) 



2^12 a;i2 



•^12 ■^12 



-2xl^,Qi/4 



'2xl^QllA 



-2's:\^,Qlli 



(3.10) 



which is the result obtained by Kovchegov in At first sight, there is a difference of the factor of 2 in the 
exponentials, which comes from effective double interaction of the same dipole. At this point we have to clarify this 
coefficient in powers of the exponentials in Eqs. (3.7)-( ^!lc| ). We found some confusion in the literature regarding 
this coefficient, it varies from 1/4 to 1/2 depending on the author and even sometimes for the same author. This 
confusion stems from the definition of the saturation scale for a quark or a gluon which differs by a factor of 2. We 



use the quark saturation scale throughout the paper having a coefficient 1/2 in the exponentials of Eqs. (3.7)-(3.10) 
The physical meaning of this is explained as follows. 



The important feature of the expression for the inclusive single gluon production Eqs. (3.7)-( 3.10 ) is that 
interaction happens only to the emitted soft gluon and a quark(antiquark) from which it was emitted. The interaction 
of a quark (antiquark) which does not emit the soft gluon cancels and can be schematically explained as follows. In 
terms of the amplitude, diagrams B{B*) and D in Fig. ^ have only two possibilities for the real soft gluon: to 
emitted either before or after the interaction time t = as shown in Fig. ^. 
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Figure 6: This illustrates the fact that in the low energy limit the interaction occurs only with the produced gluon and 
a quark (antiquark) from which it was emitted. The cancellation leads to effective interaction twice with the same (upper) 
dipole. 



The emission after r = has a minus sign w.r.t the emission before t = 0, which can be easily checked by 
considering hght-cone energy denominators (see for example [To| and also Appendix A). For simplicity, 

consider one i-channel gluon interaction with the target for each dipole (12) and (20) formed by the soft gluon 
emission. In this case we have two contributions in the transverse coordinate space for the lower dipole 



{(^ 



'ilxi 



-ilxo 



-ilX2 



)} 



-ilx2 



-ilxQ ^ 



the upper dipole has only one contribution 



72" 



-ilX2 



-ilxQ \ 



(3.11) 



(3.12) 



From Eq. ( 3.11 ) and Eq. ( 3.12 ) one can see that the two contributions for the lower dipole effectively sum into 
interaction with the upper dipole. It should be mentioned that one cannot sum Eq. ( 3.11 ) and Eq. ( 3.12| ) in the 
amplitude since they have different colour matrices, but after multiplication by conjugate amplitude one obtains the 



same colour factor which results into a factor of 2 in the exponentials in Eqs. (3.7)-(3.10) and looks like a dipole 
interacts twice with the target (or two dipoles of the same size interact once). In the BK equation (written for 
total cross section) we do not observe this, since the second diagram in Fig. ^ cancels with corresponding virtual 
gluon emission. This ad hoc explanation is by no means an exhaustive one, but rather a simple argument why one 
should expect only the soft gluon and the quark (antiquark) from which it was emitted to interact with the target. 
A rigorous calculation shows the same result, provided we assign to the interaction with target the same expression 
before and after the gluon emission. Here we demonstrated how this happens in the case of the eikonal interaction at 
low energy with no evolution included. In the next Section we show that proper inclusion of the non-linear evolution 
effects does not change this general feature of the single inclusive cross section that the "spectator quark" never 
interacts with the target in the course of the evolution and thus all the interaction can be expressed in terms of the 
adjoint amplitude Na{x) — 2N{x) — N'^{x) as was suggested by Kovchegov and Tuchin in Ref. 



4. Inclusive gluon production with evolution 

In this section we want to include the evolution into the expression for the single gluon inclusive cross section 
found in Section ||. It was suggested in Ref. Q that each exponential in Eqs. (3.7)-( ^TC| ) is to be replaced with a 



function 1 — Nq-, where Nq is the elastic scattering amplitude of an adjoint ("gluonic") dipole given by Na{'22') = 
2iV(22') — N'^{22') with A^(22') being a solution for the BK equation. This substitution is based on the classical 
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expression where only the produced gluon and a quark(antiquark) from which it was emitted do interact with the 
target and thus in the course of the evolution one can safely ignore the spectator quark (antiquark). Kovchegov and 
Tuchin |^ analyzed the next step of the evolution in the rapidity where an extra softer gluon is emitted and found 
that this form also holds to this order. This way the authors generalized the result to any order of the softer gluon 
emissions. In this Section we use the expression for the single inclusive cross section Refs. ( |3.2| - ^^ ) supplemented 
by the evolution equation for the function M{ij\ik) to show explicitly the validity of the adjoint("gluonic") dipole 
structure at any order of the evolution. 

We start with deriving the evolution equation for the function M{ij\ik), which has the meaning of the total cross 
section for j = k. It should be emphasized that M{ij\ik) is also a function of the rapidity though it is not indicated 
in our notation. We are interested in the evolution in rapidity where only softer gluons are emitted, the emission of 
the harder gluons can be easily included through the dipole density as it was shown by Kovchegov and Tuchin 
and thus is not relevant for our discussion. We want to remind that we are interested only in the situation where 
the measured gluon is the hardest gluon in the system. In our derivation we help with a set of the simple mnemonic 
rules, as follows. 



the gluon is represented by a double quark line in the large Nc limit, the triple gluon vertex can be effectively 
written as an emission of the softer gluon from the quark and antiquark components of the harder gluon. The 
emission from the antiquark component of the harder gluon has a relative minus sign w.r.t. the emission from 
the quark component. For more details, see Ref. [jl^ ; 

the virtual emission (diagrams R,C, E and F in Fig. ^ ) has relative minus sign w.r.t. to the real emissions 
(diagrams A, B and D in Fig. |^). In addition to this, the diagrams where the virtual gluon is not present at 
the interaction time t — has a factor of 1/2 due to light cone time ordered integration. For more details, see 



Ref. 1 14 



• the softer gluon can be emitted only after the harder gluon before the interaction time t — 0, and only before 
the harder gluon after the interaction time. In the case where one of the gluons is emitted before r = and 
the other one after t — 0, the light cone time ordering is arbitrary, irrespectively whether the soft or the hard 
gluon was emitted before r = 0. These light cone time ordering rules where derived by Jalilian-Marian and 
Kovchegov [l0[| using light cone energy denominators. 

We would like to comment more on the last light cone time ordering rule. This rule implies that a softer gluon 
can be emitted only from a dipole which was present at r = 0, this means the all the information about further 
evolution can be encoded in the function M{ij\ik) through specifying the coordinates of the dipoles present at r = 
in the amplitude and conjugate amplitude (and, of course, the energy variable not reflected in our notation). In 
the more general case, where the rapidity separation between two adjacent emitted gluons is not large (multi Regge 
kinematics does not hold) one has also to specify the dipoles present at t = oo. 

We are ready to write the evolution equation for the function M(12|12'), which has a meaning of the total 
scattering cross section for a dipole having quark at xi and antiquark at X2 in the amplitude, as well as, quark at Xi 
and antiquark at x'2 in the complex conjugate amplitude. We analyze the emission of a softer gluon with transverse 
coordinate X3 in the dipole "12" ("12"'). Working in the large Nc limit we represent the gluons by double quark- 
antiquark line and show them disconnected from other lines reflecting all possible couplings to quark or antiquark 
line in accordance with the notation used in In analogy with Fig. ^ consider all possible emissions of gluon "3" 
in the dipole "12" ("12'") depicted in Fig. ^ ( diagram B* and C* are not shown, but must be included as well). 

The evolution equation is found by accounting for all interaction possibilities of the system of the initial dipole 
with one extra softer gluon emission depicted in Fig. |^ (together with diagrams B* and C* not shown in Fig. 0). 
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Figure 7: All possible emissions of the softer gluon "3" in the dipole "12" ("12'"), diagrams B* and C* are not shown. 
The corresponding evolution equation is given by 
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[M(13|12') {1 - Ar(32)} + N{32)N{12')] 
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_ a^2'3 
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[Af(13|12') {1 - iV(32)} + iV(32)iV(12')] 
[M(32|12') {1 - 7V(13)} + 7V(13)7V(12')] 
[M(12|13) {1 - iV(32')} + 7V(12)7V(32')] 
[M(12|32') {1 - A^(13)} + N{l2)N{li)] 



(4.5) 



f X23, a:2'3 

1 / X2d. X2Z 



M(12|12') 
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with N being a solution to the BK equation. 
Each contribution is explained as follows 



Xi3 a:i3 

2 2 
^13 ^13 
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7V(12)7V(12') 



(4.6) 



a;2'3 a;i3 ^13 a;2'3 

"To 9 



*''2'3 ■''13 



^3 ■^2'3 



7V(12)7V(12'; 



The two terms on the r.h.s of Eq. (4.1) come from the diagrams R and R* in Fig. |^and represent the reggeization. 
Their Kernels are different and correspond to emission and absorption of the softer gluon " 3" in the amplitude 
for diagram R, and the emission and absorption of the softer gluon "3" in the conjugate amplitude for diagram 
R* . A factor of ^ appears due to time ordered integral in the light-cone time (for more details see Appendix 
of Ref. M). 



• Terms in Eq. (4.2) stand for the interaction of each dipole with the target as shown in diagram A. Their 
Kernel corresponds to the emission of the softer gluon "3" in the amplitude and its absorption in the conjugate 
amplitude. The last non-linear term in Eq. (4.2) describes the situation where the two dipoles scatter both 
elastically and inelastically off the target simultaneously. 



• For all the terms in Eq. (4.3) we have we same Kernel as for terms in Eq. ( |4.2| ) since emission and absorption 
of the softer gluon is the same. The first and the second terms reflect the situation where one of the dipoles 
rescatter both elastically and inelastically, and the other one rescatters only elastically in either amplitude 
or conjugate amplitude. These contributions are not present in Af(13|13)M(32|32') since M by the definition 
includes inelastic part and elastic scattering both in the amplitude and the scattering amplitude (in other words, 
only N'^ and not N) . The last term corresponds to the situation where the two dipole rescatter elastically, one 
in the amplitude and another in conjugate amplitude, and vice versa. 



The same analysis implies in Eq. (4.4) requiring the proper counting of all interaction patterns in diagrams 
B (first two lines) and B* (last two lines) in Fig. ^. The zero term in the asymmetric Kernel used here stands 
for situation where the dipole "13" ("12"') is not present by time r = in the conjugate amplitude and thus is 
omitted as was explained in more details in Section ||. The last two lines of Eq. (4.4) correspond to the diagram 
B* , which is the conjugate of diagram B in Fig. ^ 



• In Eq. (4.5) we account for diagram C(first two lines) and diagram C*(last two lines) using the same arguments 
as for Eq. (4.4). The only difference between the corresponding lines of Eq. (4.4) and Eq. ([l.5|) is the Kernel 
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of the dipole splitting, namely, in Eq. (4.4) the soft gluon 3 is present at r = oo, and in Eq. (4.5) is not. If 
the initial dipole "12" ("12"') has the same coordinates in the amplitude and the conjugate amplitude then the 



corresponding lines of Eq. (4.4) and Eq. (4.5) are subject to their mutual cancellation (real- virtual cancellations) 
and was shown by Chen and MuUer |14| . 



In Eq. (4^) we write the contributions of diagrams D, E and F in Fig. Here one should note that in the 
crossed dipole splittings ( i.e. fj^p^ etc.) the inelastic interaction are suppressed in the large Nq limit, this 
is the origin of the terms iV(12)iV(12'). 



The evolution equation should be supplemented by the initial condition Mq given by Eq. (B-12) 



Afo(12|12') = l + e-^22'Q'/4_e- 



(4.7) 



One can easily see that due to the equality of the quark coordinate xi in the amplitude and the conjugate 
amplitude, some terms cancel each other in the evolution equation Eqs. ( |4.lD -( pr5| ). For example, the second line in 
Eq. (4.4) cancels the second line in Eq. (4.5). Some further cancellations happen because of the optical theorem. 
The optical theorem states that the total cross section equals twice the imaginary part of the scattering amplitude 
M(12|12) 



2iV(12). This implies that the last term in Eq. (O) cancels the first term in Eq. (4.3) 



The simplified form of the evolution equation Eqs. (4.1)-(4.5) is given by 
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7V(13)iV(32') - 7V(32)7V(13)} 



(4.9) 
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[M(13|12') {1 - 7V(32)} + A^(32)iV(12')] 



(4.10) 



^_^) Af(12|12') 

^23 -''2'3- 



(4.11) 



To see the consistency with the previous studies one can take 2 = 2' and use the optical theorem to recover the 



BK equation. In this case Eq. (4.10) and Eq. (4.11) vanish. Kernels in Eq. (4.8) and Eq. (4.9) become equal and we 
end up with with the evolution equation 



d{2N{l2)) ^ 
dy 2tt 





X2i\ 


/ (fxs, { ^ 




/ Wis. 


^23 / 



{27V(13) + 2iV(32) - 27V(13)^(32)} 



with the initial condition directly obtained from (see Eq. ( B-12 )) 

Mo(12|12') l + e"^22 



/4 _ g-^^2Q.o/4 _ g- 



,/4 



(4.12) 



(4.13) 
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for 2 = 2' as follows Mo(12|12) = 2 (l 
with the correct initial condition. 



L/4 



2A''o(12). Thus the evolution equation is the BK equation 



( [4.111) . At this 



Next step is to study general properties and to find a solution to the evolution equation Eqs. (4. 
point we make a digression and go back to the classical expression for the single inclusive production cross section 
Eqs. ( |3.2D -( |3^ . As one can easily see, the cancellation of the emissions with the "spectator" quark (antiquark) found 
byKovchegov g in the form of Eqs. (^-(|3.10D can happen only if A'fo(12|10) = No{12) + No{10) - No{20) , where in 
the last term 20 — (10) — (12). The initial condition Eq. ( [4.13 ) indeed satisfies this condition. If one wants to include 
the evolution effects, the function Mq in Eqs. (3.2)-( p^ ) should be replaced by the solution to the evolution equation 
Eqs. (4.8)-(4.11) with the initial condition given by Eq. ( 4.13| ). This suggests that M(12|12') should also have this 



property to result into the cancellation of the emissions and the interactions with the "spectator" quark (antiquark) 
line. The straightforward substitution of 



M(12|12') = iV(12) + N{12') - N{22') 



into Eqs. ([4.8[)-( [4.1l|) reveals 

a(Af(12) + iV(12') - jV(22')) 
dy 
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2^13 



X23 



Xl3 
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(Ar(12) +iV(12') -iV(22')) (4.15) 
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{2iV(13) + (iV(32) + iV(32') - N{22')) - Ar(13)iV(32') - Ar(32)iV(13)} (4.16) 



xi3 _ £23 \ ^£23 _ x^\ [(•jv(13) + iV(12') - iV(32')) {1 - N{32)} + N{32)N{12')] (4.17) 

^13 ^23/ \^23 ^2'3 / 



2'3 -^23 



X2'3 X23 \ I Xi3 X2'3 



] [(iV(12) + iV(13) - iV(32)) {1 - Ar(32')} + iV(12)Ar(32')] 



''2'3 

1 f X23 X 2'3 
2 

2' 3 



2 V a;23 ^'o' 



(7V(12) +iV(12') -iV(22')) 

which is just a linear combination of three BK equations for initial dipoles with coordinates 12, 12' and 22' 



(4.18) 



d{N{12) + N{12') - N{22')) _ a, 
dy 27r 




d'xzli^f-^] {iV(13) + iV(32)- Ar(12) -7V(13)iV(32)} (4.19) 



2^13 



2 

^13 



2 

{A^(13) + iV(32') - iV(12') - iV(13)Ar(32')} (4.20) 



{iV(32) + 7V(32') - 7V(22') - iV(32)7V(32')} j. (4.21) 

■^23 ■^2'3/ J 

This way we show explicitly that indeed the cancellations claimed by Kovchegov and Tuchin |^ hold to any order 
of evolution, allowing for the simple description of the single inclusive cross section in terms of the combination 
2N {x) — N"^ [x) , which can be interpreted as the scattering amplitude Ng{x) of an adjoint ("gluonic") dipole of size x. 

In the next Section we examine the Abramovsky-Gribov-Kancheli cutting rules in pQCD using the formalism we 
developed for the calculation of the single production cross section. 
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5. Abramovsky-Gribov-Kancheli cutting rules 



More than 30 years ago Abramovsky,Gribov and Kancheli |^ analyzed the relative contributions of processes with 
different multiplicity of the produced particles to the total cross section of the multiple Pomeron exchange in the 
framework of the Reggeon Field Theory. They found that the contributions of different unitarity cuts, corresponding 
to different particle multiplicities are proportional to each other (AGK cutting rules), and their relative coefficients 
lead to some peculiar cancellations, called AGK cancellations for the single jet production. The most interesting case 
is the two Pomeron exchange, where the AGK cancellation can be explained as follows. The total cross section is 
given by the sum of all possible unitarity cuts. The Pomeron is non-local object, so it can be also cut enlarging the 
number of cut possibilities. This way one obtains two distinct objects - cut and uncut Pomerons, which are to be 
treated separately as far as multiplicity is concerned. The simplest example for this is the two Pomeron exchange 
depicted in Fig. ^. As one can see we have here three different processes: a) double cut, where the two Pomerons are 




total double cut single cut diffractive cut 

cross section 



Figure 8: The total cross section of two Pomeron exchange in terms of contributions from different s-channel unitarity cuts. 
The double cut diagram describes the process with the multiplicity two times larger the multiplicity for one Pomeron exchange 
(see Fig. |l|), the single cut diagrams leads to a process with the same multiplicity as one Pomeron exchange (see Fig. |l]) and 
the diffractive cut diagram gives the expression for a diffractive production with small multiplicity. 

cut. This brings two times multiplicity of one cut Pomeron; b) single cut, one Pomeron appears cut and the other 
one is uncut. The multiplicity of this process equals the multiplicity of one cut Pomeron; c) diffractive cut, here 
two Pomerons are uncut and no particles are produced. The multiplicity of this process is very small. The relative 
weights of three processes are calculated as follows. There only one possibility to cut two Pomerons at the same 
time, each cut Pomeron is twice the uncut Pomeron exchange ( due to the unitarity condition f = 2P ), thus we 
obtain AP^ . For the single cut we have two possibilities: to cut one of two Pomerons, this must be multiplied by 
the a factor of 2 for putting the uncut Pomeron either in the amplitude or the conjugate amplitude. The resulting 
combinatorial coefficient of 4 is to be multiplied by a factor 2 for the normalization of the cut Pomeron {f = 2P) as 
well as a minus sign originating from the fact that uncut Pomeron is the imaginary part of the scattering amplitude. 
Summarizing all this, one gets —8P^ from the single cut. Finally, in the diffractive cut the combinatorial coefficient 
is 2 for interchanging the two Pomerons in the amplitude and the conjugate amplitude. All three cuts add up to 
4P2 - 8P2 + 2P2 ^ _2p2 ^ ^(2)^ as expected. 

This set of the AGK coefficients results into cancellation of the single jet production as shown in Fig. ^ 

The jet can be produced from either of two cut Pomerons in the double cut, and only from the cut Pomeron 
in the single cut, resulting into the cancellation of the single jet production in the two Pomeron exchange. In this 
Section we treat the contributions from the different unitarity cuts in pQCD in the framework of the color dipole 
model. The goal of this study is to see how the AGK cancellations happen in the color dipole model and check in 
this way the validity of the AGK cutting rules in pQCD. 
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Figure 9: The cancellation of the single jet production for two Pomeron exchange. 



In the original derivation of the AGK cutting rules [g| some basic assumptions were made, namely, 

• Pomeron is a proper asymptotics in the high energy scattering; 

• there are no "half-cut" Pomerons, the contributions where a Pomeron is not cut entirely are suppressed expo- 
nentially in the invariant mass; 

• there is no contribution to the multiplicity from the cut of the vertices; 

• all the vertices for various cuts are the same and real. 

In the AGK paper |^ all these assumptions have been proved in the superconverged field theories, since they guarantee 
that all scattering amplitudes decrease in the region of large masses. As has been shown in Ref. in the leading 
log{l/x) approximation of perturbative QCD, the amplitudes fall in the region of high mass fast enough to hope that 
AGK cutting rules are correct. However, in addition to the assumptions that discussed above, it was also assumed 
in Ref. Q that the multiparticle production is described by the same set of diagram as the total cross section. This 
assumption was doubtful even in a superconverged theory, since it was found that there is an example of the diagrams 
that does not contribute to the total cross section, but could lead to the multiparticle production ( the AFS-type 
of diagrams [p6| ). Therefore, our main goal is to find what types of diagrams indeed contribute to the total and 
inelastic cross sections and to calculate the relation between different production processes. 

Our strategy is as follows. We use the formalism developed in the previous Sections for deriving cutting rules in 
pQCD, then we analyze the single inclusive production of the two Pomeron exchange and check how the emission of 
the gluon changes the cutting rules. The process we discuss is the dipole scattering on the target via two Pomeron 
exchange: a) two Pomerons couple to the dipole; b) one Pomeron couples to the dipole and then splits into two 
Pomerons. The first case trivially satisfies AGK cutting rules and thus is of no interest. We focus on the second 
case depicted in Fig. At this stage we do not consider any gluon production and want to analyze only the two 
Pomeron exchange in terms of the contributions of different cuts. 

Schematically, this process can be written as 



where P{Y — Yq) is the Pomeron amplitude and Pap is the triple Pomeron vertex. In the case of the virtual photon 
scattering on the target this is given by 




(5.1) 
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Figure 10: We check the vahdity the AGK cutting rules for the process of one to two Pomeron sphtting. 



cr = J d^xid^x^dz d^xid^xo l^p'^ '^'^'^ {xiq, z)\'^n{xi, Xq\xi, xo;Y — Yo)a2p{xi, xq (5-2) 

where ip'' ~*'^(a;jQ, z) is the wave function of the virtual photon being spht into qq pair, n{xi, Xq\xi, xq; Y — Yq) is the 
dipole density, which is the solution to the linear evolution equation. The triple Pomeron vertex and the interaction 
with the target is encoded into cr|p(a;i, xq). Using the formalism developed in the previous Sections we can readily 
write an expression for cr|p(a;i,Xo) as 



a|f,(xi,a;o) = ^ ^ dy f d^X2 - ^) (M,™(12|12)Mi,„(20|20) + i (M,,„(12|12))' 

+ i (Mh„(20|20))' - i (M,,„(10|10))' - {AfH„(12|12) + Mi,„(20|20)}{27VBFKL(12) + 2iVBFKL(20)} (5.3) 
+Mi,„(10|10)2iVBFKL(10) + 2iVBFKL(12)iVBFKL(20) + (A^bfkl(12))' + (iVBFKL(20))' - (7Vbfkl(10))' } 

where A'bfkl and Alun are solutions to the BFKL equation and some linear evolution equation as explained below. 

In Eq. (5.3) the triple Pomeron vertex is represented by the initial dipole "10" which splits into two new dipoles 
" 12" and "02" , each of them can interact with the target elastically by the BFKL amplitude Nbfkl and inelastically 
by the inelastic cross section Mun . If the same dipole interacts twice inelastically we put a factor of ^ for the identity 
of the interactions. We have to take into account all possible splittings of the initial dipole as shown in Fig. since 
most of the splittings are canceled by the real- virtual cancellations and we are left with diagrams A, R and R* . The 
inelastic cross section Mun is the solution to the linear evolution equation, obtained from the non-linear evolution 
equation Eqs. (4.1)-(4.6) retaining only the its linear part. In the contrary to the the non-linear case we do not include 
multiple rescattering in the initial condition for Mun, which now consists of only two t-channel gluon exchange. In 
this case the total cross section equals the inelastic cross section and Mun is associated with a cut Pomeron. For 
the non- linear case, where the multiple rescattering is taken into account the solution to Eqs. (4T)-(4^) includes 
both cut and uncut Pomeron exchanges, which can be extracted by expanding the total cross section in terms of the 
Pomerons as we show below. As we have already mentioned the linear evolution equation for Mun is readily obtained 
from Eqs. (|4.l| )- (|4.6|) neglecting the non-linear terms as 

^^^^4^ = 1^ I d^^. l-l - M.„(12|12') ^)^/.„(12|12') (5.4) 
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with the initial condition given by 

Mh„o(12|12') 



2^23 2^2'3 

2 2 
■^23 ^2'3 



^^)Af,,„(12|12') 

2^13 ^13 / 
z V J-23 ■''23 2^3 J.]^3 / 



2^2'3 a^2'3 

2 2 
2^2'3 ^2'3 



5^^)Af.„(12|12')| 

■''13 2-13 / J 



)(e 



ikxi 



fc4 



(5.8) 



-/F(fc2 



(5.9) 



where IF(k'^) is the impact factor of the target. 

The terms on r.h.s of Eq. (5.4) are the contributions from the diagrams R and R* in Fig. §, Eq. (U) comes from 
diagram A, Eq. ( |5.6|) account for diagrams B and B*, Eq. (5/7) stands for diagrams C and C*, and, finally, Eq. ( |3.6[ ) 
describes diagrams Z?, £' and F. A more detailed explanation is presented in Section ^. 

Eqs. (^.4|)-(^^) can be written in a compact form as 
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For X2 = X21 all terms in Eqs. ( 5.12| )-( 5.12 ) cancel out and the linear evolution equation reduces to 
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;Mi,„(12|12) + Mi™(13|13) + AfH„(32|32) 



(5.13) 



which is BFKL equation because of the unitarity condition Afii„(12|12) = 2Niin{12). As in the case of the non-linear 
equation it is easy to see from Eqs. (5.10)-(5.12) that its solution A//i„(12|12') can be expressed through solution to 
BFKL equation as 



Afi,„(12|12') = NBFKLil2) + NbfKl{12') - NBFKLi22') 



(5.14) 



This is a generalized form of the optical theorem, which relates the total cross section Af/i„(12|12') of scattering of 
the dipole "12" ("12"') to the elastic scattering amplitude of the dipoles "12", "12'" and "12 - 12'". 

We are interested in the contributions to the total cross section, which have different multiplicities as shown in 
Fig-O 






Figure 11: The contributions to the total cross section having different multiplicities. 



We go back to the expression for the triple Pomeron vertex with interaction Eq. (5.3) and read out the contri- 
butions of the different cuts as follows. 



Double Cut 
Single Cut 
Diffractive Cut 



i (AfH„(12|12) + Munm20)f - iAi-,,„(10|10) 
(Af,,„(12|12) + A//,,„(20|20)) (27Vbfkl(12) + 2iVBFKL(20)) + A.fH„(10|10)2ArBp,L(10) (5.15) 

1 (iVBPKL(12) + iVBFKL(02))' - |iVB'FKL(10) 



One can easily see using the unitarity relation Eq. (5.14) Miin{ij\ij) — 2iVBFKL(ij) both the double and the single cut 
contributions are proportional to the diffractive term as expected from the AGK cutting rules. 
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Double Cut 
Single Cut 
DifFractive Cut 



2 (A^BFKl(12) + iVBFKL(02))' 

-4(7Vbfkl(12) + 7Vbfkl(02))' 



2A^Ifkl(10) 

MiVlFKL(lO) 



(5.16) 



(7Vbfkl(12) + A^bfkl(02))' - Nl^^{10) 



It is instructive to demonstrate how terms in Eq. ( 5.16| ) add up to the total cross section. The easiest way to see this 
is to write down the Glauber expression for the total cross section in Eq. ( |5.3|) as follows 
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where we used the notation of Appendix H and absorbed for simplicity r(6; i?^) in the definition of a . The 



first term in the brackets in Eq. ( 5.17 ) is the total cross section of the scattering of qqg system after the soft gluon 
emission (dipole splitting) happened, this term comes from diagram A if Fig. |. The second terms in Eq. ( p.l7[ ) 
corresponds to the situation where the gluon is both emitted and absorbed before the scattering takes place (see 
diagrams R and R* in Fig. |^). It is easy to see that the expansion of the expression in the brackets in Eq. (5.17) to 
the second order in a^^ with the help of — 2N^^ gives the sum of terms in Eq. ( ^.16 ). 



The coefficients in Eq. ( 5.16 ) show that we fully reproduce the AGK cutting in this case. However, the situation 
changes when one emits a gluon from the triple Pomeron vertex. In this case the real- virtual cancellations does not 
happen and one should expect non-trivial contributions from different cuts as discussed in the next Section. 



6. AGK rules for a production from vertex 



In deriving the triple Pomeron vertex Eq. (5.2) we used the real-virtual cancellations, where the real gluon emissions 
were cancelled by the virtual emissions. This does not happen for the real gluon production (see Fig. |l^) and we 
have to go back to extract only the real emission terms. 




Figure 12: The emission from the vertex in one to two Pomerons splitting. 



As in the case of no gluon production, this process can be written as 

da 



dYo d'^k 



P{Y - Fo) ® r3p(fc, Fo) ® P\Yo - 0) 



(6.1) 
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where P{Y — Yq) is the Pomeron ampHtude and T^p is the triple Pomeron vertex with one gluon emission. In the 
case of the virtual photon scattering on the target this is given by 



da 



dVo d^k 



= / d^xid^x^dz d^xid^XQ l?/'^ ^''''(a;j5, z)pr7,ii„(a;];, XqI^i, xo; F - Fo)- 



d(7fp(xi,xo) 



dYod^k 



(6.2) 



where ip'^ ^''9(a;jg, z) is the wave function of the virtual photon being split into qq pair, niin{xi, Xq\xi, xq, Y — Yq) is 
the dipole density, which is the solution to the linear evolution equation. 



All the information about the contributions of the different unitarity cuts is encoded in the differential cross 
section — gfos^" ■ We can readily write it following the strategy we used before in deriving the single inclusive cross 
section Eqs. (3. 2)- (3. 6) in Section ||. 



+Mz„(12|12')M,,„(20|2'0) - (Mh„(12|12') + AfH„(20|2'0))(iVBFKL(12) + iVBFKL(12') + AfBFKL(20) + iVBFKL(2'0)) 

+iVBFKL(12)^BFKL(2'0) + iVBFKL(20) ArBFKL(12') + iVBFKL(12)^BFKL(12') + A^BFKl(20)7Vbfkl(2'0)} 



^ - ^ ) ( - ^ ) {-M;m(12|10) (7Vbfkl(20) + 7Vbfkl(12) + iVBFKL(lO)) 

^12 ^02 / V ^02' / 

+ iA/,2„(12|10) + iVBFKL(12)iVBFKL(10) + AfBFKL(20)7VBFKL(10)} 

a;i2 ^02 \ (^^_q\ {_Mi,„(20|10) (AfBFKL(12) + AfBFKL(20) + A^bfkl(IO)) 



^19 



"12' 

1 9 
- f2 



+ -M,i„(20|10) + iVBFKL(12)iVBFKL(10) + AfBFKL(20)7VBFKL(10)} 
- ^) - ^) {-Mu„il0\l2') (7VbFKl(2'0) + iVBFKL(12') + iVBFKL(lO)) 



X, 



02 



''12 





2^02' A 








^02' / 








X02'\ 


) 


^02' / 



¥ - ) ( ^ - ^ ) {-Mi»n(10|2'0) (iVBFKL(12') + iVBFKL(2'0) + iVBFKL(lO)) 



+ ^^L(10|2'0) + iVBFKL(10)AfBFKL(2'0) + iVBFKL(10)iVBFKL(12')} 



Xi2 Xi2' , X02 X02' \ f 1 , ,r2 



•^12 -^12' ■^02 ■''02' '' ^ ^ 
Xi2 X02' , X02 Xi2' 



Xl2 Xq2' Xq2 Xi2' 



-M/„(10|10) - M,„(10|10)2iVBFKL(10) + iVB^FKL(lO) 

^^Ifkl(IO) 
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We read out the double cut contribution from Eqs. (^.3|)-(6.4) 
dYod^k 27T (27r)2 y « ^2 « X2'e 

=^^12:1 {(m2^^(i2|i2') + mL(20|2'0))2 - ML(10|2'0) - ML(20|10) + AC(10|10)} 



3^02 Xo2' 1 
•^02 •*'02' 

{(ML(12|12') + ML(20|2'0))2 - ML(10|2'0) - ML(12|10)} 

■^12 -^02' ^ 



i- {(M2„(12|12') + ML(20|2'0))2 - ML(10|12') - AC(12|10) + ML(10|10)} 



+ ¥T J {(ML(12|12') + ML(20|2'0))2 - ML(10|12') - ML(20|10)} 
It is a tedious, but a straightforward calculation to find the single cut contribution in the form of 

d<^2Pis^ngle)(^l'^0) _ as 1 



dYod^k 2tt (27r) 



d^X2 d'^X2'e~'''^'-'^-'-'^'^ (6.5) 



{- (Mh„(12|12') + A/,™(20|2'0)) (iVBFKL(12) + A^bfkl(20) + iVBFKL(12') + N^vkl{2'Q)) 

+Mh„(20|10) (7Vbfkl(10) + iVBFKL(12) + AfBFKL(20)) + A/h„(10|2'0) (TYbfrlIIO) + iVBFKL(12') + N^fkl{2'Q)) 
-2A^h„(10|10)7Vbfkl(10)} 

2^12 X02' 



{- (A//,,„(12|12') + AfH„(20|2'0)) (iVBFKL(12) + A^bfkl(20) + A^bfkl(12') + iVBFKL(2'0)) 



12 -^02' 

+AfH„(12|10) (7Vbfkl(10) + iVBFKL(12) + AfBFKL(20)) + A^f,„(10|2'0) (TVbfklIIO) + A^bfkl(12') + AfBFKL(2'0))} 

_^X12^ {- (A^,„,(12|12') + A'/h„(20|2'0)) (iVBFKL(12) + iVBFKL(20) + A^BFKl(12') + AfBFKL(2'0)) 
2^02 2^12' 

+ AfH„(20|10) (iVBFKL(lO) + iVBFKL(12) + AfBFKL(20)) + A/h„(10|12') (7VbfKl(10) + iVBFKL(12') + AfBFKL(2'0))} 

+ {_ (A^;,„(12|12') + A./h„(20|2'0)) (iVBFKL(12) + iVBFKL(20) + A^BFKl(12') + AfBFKL(2'0)) 

^02 ^02' 

+ AfH„(12|10) (iVBFKL(lO) + iVBFKL(12) + AfBFKL(20)) + Mh„(10|12') (7Vbpkl(10) + iVBFKL(12') + AfBFKL(2'0)) 

-2A^h„(10|10)A^bfkl(10)} 
Finally, the diffractive cut gives 

rfsVin(^i>^o)_a. 1 /^2^^^2^^,^-.fe(..-.„)f^„^]f^_:^£Hl) (6.6) 



dYod'^k 2-K (27r)2 j \X\2 Xq2 

{iVBFKL(lO) - iVBFKL(12) - A^bfkl(20)} {A^bfkl(IO) - A^bfkl(12') - A^bfkl(2'0)} 



The contributions from the different cuts can be further simplified using the unitarity condition Eq. (5.14), 
but the expressions become cumbersome and we want to consider only a special case of ^^^^ dipole splitting for 

^12 ^12' 

illustrating the general conclusion regarding the AGK violation. The cuts are: 
• double cut 

^(4iVBFKL(10) - (A^bfkl(IO) - iVBFKL(12) + N^YYc^mf " (A^bfkl(IO) - A^bfkl(12') + A^bfkl(2'0))2 (6.7) 

2a 



+ (iVBFKL(12) + A^bfkl(12') + iVBFKL(20) - 2A^bfkl(22') + A^bfkl(2'0))2) 



• single cut 



- 2(iV|pKL(10) + A^Ifkl(12) + (A^bfkl(12') + A^bfkl(20))(7Vbpkl(12') - A^bfkl(22')) (6.8) 

+ (7VbFKl(12') + A^BFKl(20) - iVBFKL(22'))A^BFKL(2'0) - iVBFKL(10)(^BFKL(20) + 7Vbfkl(2'0)) 
+A^BFKl(12)(A^BFKl(12') + iVBFKL(20) - A^BFKl(22') + iVBFKL(2'0)) 

• diffractive cut 



{^bfkl(IO) - A^bfkl(12) - A^bfkl(20)} {iVBPKL(lO) - iVBFKL(12') - A^bfkl(2'0)} (6.9) 



It is clearly seen from Eqs. ( |6.7[ )-(p.9|) that already in the single inclusive cross section the original AGK cutting 
rules are violated in the sense that the contributions from the different cuts are not proportional to each other. This 



should be compared to the cuts contributing to the total cross section given in Eq. (5.16), which are proportional to 
the total cross section with coefficients 2 — 4 ~ 1. 

It should be emphasized that, in contrary to the original AGK cutting rules, non of the cuts is expressed through 
the total cross section as we show in Appendix |^. 



In the production of another gluon from one of the lower Pomerons (see Fig. 13 ) one would expect a cancellation 



Figure 13: The emission from the vertex in one to two Pomerons splitting, 
naively applying AGK cutting rules as shown Fig. M. This does not happen if the gluon production is preceded by the 



emission from the vertex, because in this case the double cut Eq. (6.7) is not proportional to the single cut Eq. (6 

Another important remark is to be made, in the real gluon production from the vertex we took into account only 
diagrams A, B, B* and D in Fig. g As we showed in Section ^ the cuts of diagram A fully satisfied the AGK rules 
for the total cross section, and we expect the AGK violation to originate only from diagram B, B* and D in Fig. |[ 
This is not the case, because of the fact that the produced gluon transverse coordinate in the amplitude X2 differs 
from its coordinate X2' in the conjugate amplitude. The cuts of diagram A in Fig. for the dipole splitting ^^iiai 

are given by 

• double cut of diagram A 

\ (AfBFKL(12) + AfBFKL(12') + AfBFKL(20) + AfBFKL(2'0) - 2AfBFKL (22') ) ^ (6.10) 
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• single cut of diagram A 

- (AfBFKL(12) + iVBFKL(12') + AfBFKL(20) + ArBFKL(2'0)) (AfBFKL(12) + iVBFKL(12') + JVbFKL (20) + iVBFKL(2'0) - 2AfBFKL (22' )) (6.11) 



• difFractive cut of diagram A 

+ (AfBFKL(12) + iVBFKL(20)) (iVBFKL(12') + AfBFKL(2'0)) 



(6.12) 



The double Eq. ( 6.10| ), the single Eq. (6.11) and the diffractive Eq. (6.12) cuts are proportional to each other only 

for X2 = X2' ■ 



Eq. (6.7) - Eq. ( |6.9[ ) does not allow us to calculate the double inclusive cross section but we can easily to calculate 
the following observable / (Pki(Pk2 f^^ (iy2rf^o'/(i^fci(i^fe2cJyi(iy2 which gives the multiplicity of particles with rapidity 
less than yi which accompany the produced jet at rapidity yi integrated over possible transverse momentum. From 
Eq. (3/7) - Eq. (6^) one can conclude that 



(fki(fk2 



dy2 



<Pki(Pk2dyidy2 



(6.13) 



as 1 
2tt 2tt 



(fx2 



■47Vbfkl(02) (7Vbfkl(01) - A^bfkl(02) + iVBFKL(12)) (|n(yi - 0)|)j 



where {\n{yi — 0)|)p is the average multiplicity in the rapidity window yi — 0. 



One can see that cross section of Eq. (3.13) is not equal to zero and therefore, this equations confirms the main 
resuh of Ref. 0. 

In the next Section we show processes of different multiplicity can be extracted directly from the Glauber 
expression for the single inclusive cross section. 



7. AGK rules in Glauber Formalism 



In this Section we use the Glauber expression for the single inclusive cross section including multiple rescatterings 
to find the contributions of different multiplicities to the total cross section. In the classical expression for the 



single inclusive cross section Eqs. (p.7D-( p.lO| ) we included both elastic and inelastic rescatterings on a large target 
consisting of an infinite number of nucleons. In this Section we consider only two nucleons and show that in the 
leading order this coincides with the expression for two Pomeron exchange in Section |^. In order to do this we need to 



separate elastic and inelastic contributions in the Glauber expression Eqs. (3.7)-(3.10), since in its derivation we use 
cancellations of these two as follows. In Section || we started by defining cross section MQ{ij\ik), which includes both 
elastic and inelastic parts and was found by Glauber summation in Appendix In passing from Eqs. (3.2)-( ^!6| ) 
to Eqs. (^-(p^ we used Eq. ^ 



obtained from Eq. ( [B-5| ) with the help of Eq. ( B-6 ). In this step we canceled 
some of the inelastic terms with the elastic ones. Our goal is to reconstruct the classical expression for the single 
inclusive cross section before we make use of these cancellations. To do this we just substitute the expression for the 
cross section Mo(jj|z/c) given by Eq. (B-5) in Eqs. (^)-(p^) 
After some lengthy calculations we obtain 



1 



-i(<TSA(12)+„ 



M20)) 



1 {^^ 



d?kdy 
^(12')+<T-S'*(2'0)) 



27r (27r)' 







( ^12' 


_ ^02^\ 












^02 / 


U?2' 


^02' / 



-i((T^-*(12)+CT^-*{20) + CT-^'^(12')+o--S-*(2'0)) + CTi„(12|12') + CTi„(20|2'0) 



(7.1) 
(7.2) 



/ X12 


3:02 \ 








V^12 


^02 / 





*(12)H 



'(20)) 



HlO) 



-i(<T^^-*(12) + o 



M20)4 



,(12|10)A 



(7.3) 



- 23 - 



X02 
■''02 



1^12' 
^•12' 



-i(cT«^(12)+„ 



1(20)) 



(l2)^ 



M20)H 



HlO))4 



,(20|10) 1 



- 0- 



X02 
'Z2~ 
■^02 



^12' 
■^12' 



^02' 
■^02' 



^10) 



-i(<TS'*(12')+a-B'^(2'0)) _^ g- 1 (<tS-*(12') + <tS-4(2'0)+o 



'(10))H 



,(10112')'| 



X12 



^12' 
^2 



3^02' 
-.2 



^12 Xi2i 
2 

12' 



■^12 



^02 3^02' 
''02' 



"^02 ■ 



1 - 2 e" 



'(10) 



'(10) , 



'{l2)^ 



'(10)4 



'(20)) 



,(10|10) ^ 



-i((T^'*(12') + (T-^^(2'0)+(T-^'^{10)) + CTi„{1012'0)^ 



3^12 Xq2' Xi2 
9 9 9 

02' 02 -^12 



r2 
42 



i)(l-e-^'^"-^(l«))' } (7.4) 



where, for simplicity, we absorbed the profile function T(6; Ra) in the definition of the scattering cross sections 
and ain ■ The same expression can be directly applying Glauber formula as shown in Appendix |^. 

Now, we expand the exponentials Eqs. (7.1)-(7.4) in powers of ct™ to the second order to find the contributions 
that correspond to the double cut of two Pomeron exchange found in Eq. ( |6.4D . After some algebra we obtain the 
same structure we have in Eq. (3.4). Expansion in powers of (a^^)"^ reproduces the diffractive cut in Eq. (6.7), and. 



finally, the expansion in powers of (7'"cr corresponds to the single cut in Eq. (6.5). It is important to note, that 



the classical expression in Eqs. {7A_)-{7A) does not account for the evolution effects, and thus should not reproduce 
Eqs. (7.1)- (7.4) in more general case, when the evolution is switched on. The way to do this is to go back to 
Eqs. (3.2)-(3.6) and substitute solutions to the evolution equations M{ij\ik) and N{ij) in place of MQ{ij\ik) and 
No{ij) as described in Section ^. Then, expand to the desired order in am and as A to obtain contributions of 
different multiplicities. However, because of the fact that M(ij\ik) and N{ij) are solutions to non-linear equations 
with "Pomeron splittings" where "cut Pomeron" can split to two "uncut Pomerons" at some rapidity, fixing the 
interaction type {am or as a) by expanding the solution bring little or even no information about the multiplicity 
of the produced particles in the rapidity window between the target and the projectile. Thus, the proper way for 
studying the multiplicity distribution is to attach Pomerons (cut or uncut) to a triple Pomeron vertex at some given 



rapidity, as we did deriving Eqs. (6.3)-(S.4) 



8. Conclusions 

The central result of this paper is the proof of the validity of the Abramovsky-Gribov-Kancheli cutting rules in pQCD 
for the total cross section, and the demonstration of their violation for a particle production from the vertex. This 
violations happens already for one gluon production from the triple Pomeron vertex as shown by explicit calculation 
of contributions from different cuts to the inclusive cross section of two Pomeron exchange. The general conclusion 
can be formulated as follows: AGK rules valid for 

• total cross section; 

• production of any number of gluons from Pomerons either before or after triple Pomeron vertex. 
AGK rules violated for 
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• one and more gluon productions from the triple Pomeron vertex. 

This result allows to build a generating functional that will incorporate all multiparticle production processes in the 
spirit of Ref . , and find the equation for diffractive dissociation of Ref . . 

As a by product of our analysis we derived a generalized form of the Balitsky-Kovchegov evolution equation for 
a non-diagonal cross section of the dipole scattering for a dipole having different coordinates in the amplitude and 
the conjugate amplitude. This generalised BK equation is easily solved by a linear combination of the solutions to 
the BK equation. This particular form of the solution preserves the adjoint (dipole) structure of the single inclusive 
cross section found by Kovchegov and Tuchin to any order in the evolution. 

The generalized BK equation found in the present paper is extremely useful for a systematic treatment of 
multigluon productions, multiparticle correlations including jet production in diffractive dissociation and multiplicity 
distributions in more general cases. 

Acknowledgements 

We are grateful to Jochen Bartels, Sergey Bondarenko, Errol Gotsman, Yura Kovchegov, Lev Lipatov, Cyrille Mar- 
quet, Uri Maor and Kirill Tuchin for fruitful discussions on the subject. Our special thanks go to Yura Kovchegov 
who found a mistake in the previous version of this paper which led us to an incorrect violation of Kovchegov- Tuchin 
formula for single inclusive cross section. One of us (A. P.) thanks Nestor Armesto and Carlos Pajares for their hospi- 
tality at University of Santiago de Compostela during the work of the paper. A. P. is also thankful to the organizers 
of the Les Houches School on Hadronic collisions 2008, in particular, to Yuri Dokshitzer, Francois Gelis and Edmond 
lancu for their interest in the results of this study presented at the school. 

This research was supported in part by the Israel Science Foundation, founded by the Israeli Academy of Science 
and Humanities, by BSF grant # 20004019 and by a grant from Israel Ministry of Science, Culture and Sport and 
the Foundation for Basic Research of the Russian Federation. 



- 25 - 



Appendices 



A. Light Cone Perturbative QCD at high energy (simple diagrams) 



In this Appendix we calculate simplest diagrams in the Light Cone Perturbation Theory (LCPT) applied to QCD, in 
particular, we show how the minus signs leading to the real-virtual cancellations appear in this framework. The rules 
for LCPT in QCD were formulated Ref. LCPT is very useful for high energy scattering since it represents an 
intuitive space-time picture of the scattering process. As a simple example of LCPT calculations consider a diagram 
where a soft gluon k is emitted from the heavy quark, and then the heavy quark interacts with the target as shown 
in Fig. The interaction with the target is mediated by the gluon m and the intermediate states are denoted by 
the vertical dotted lines with light-cone energy denominators Di and D2 are calculated from 



(A-l) 



where J^i -P- ^^'^ J2fP- ^"^^ sums of the light-cone energies of the intermediate and the initial (or equivalently 



final) states, correspondingly. The denominators Di and D2 are readily found using Eq. (|A-l|) 



Di = {p- 4- fc_ + P_ - p_ - P_ 



P+ 



= k- 



(A-2) 



D2 — {p ^ k — m)- 



^ ^ + ^ + tA^^^^=m^ (A-3) 



(p — k — mY 



where we used the fact that p+ 3> ^ m+ in accordance with the Rcgge kinematics. 




p-k-m 



P+m 



Figure 14: Emission of gluon in Light Cone Perturbative QCD . 



The expression for the diagram depicted in Fig. is given by 



u{p-k-m) u u{p-k) 1 
=9t 7 



n{p-k) g^a^.^ 
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After some simplifications due to the normalization properties of the spinors in the kinematic regime p+ ^ ^ m+ 
we obtain 



(A-5) 



where we used the fact that e • fc = = e+fc_ + e_fc+ — • fc^ ~ e_fc_|_ — • fc^. 

Next, we calculate the so-called late emission diagram where the emission of the gluon k occurs after the inter- 
action as shown in Fig. |l^ 















target 







p-k-m 



P+m 



Figure 15: Late emission diagram in LCPT. This diagram is opposite in sign w.r.t the early emission in Fig. |l4| 
As in the case of Fig. |l| we write the light-cone energy denominators 

D, = ip-m)^+ m_ + P_ - p_ - P_ = {^"""fi +!!A + ^^!!±=rn^ (A-6) 

[p — m)+ 771+ p+ 777+ 



D2 = {P ^ m)- + 777- + (P + ?7T,)_ — (p — k — m)- — (P + 77l)_ — k- — fc_ 



(A-7) 



In Eq. (A-7) we used the equivalence of the initial and the final light-cone energies. Note that due to the presence 
of the emitted gluon k in the final state D2 generates a minus sign. The rest of the calculation is straightforward 



and in full analogy with Eq. (A-4). The final answer for the late emission diagram in Fig. [13 is given by 



AgH" (g) t" ^2gt'' 



k- 



4gh' < 



kl 



(A-8) 



and differs from Eq. (|A-5| ) only by a minus sign coming from the light-cone energy denominator Eq. ( A-7 ). 

Here we presented a simplified case only one quark scattered off the target. In general, any late emission diagram 
will generate a minus sign w.r.t early emissions. 



B. Calculation of Mo{il\kl) 

In this Appendix we derive the explicit expression for MQ(il\kl) for dipole- nucleus scattering in the Glauber approach 
assuming the Born Approximation of perturbative QCD (the exchange of two Coulomb-like gluons) for the dipole- 
nucleon interaction. It is well known that the Glauber formula for the imaginary part of the scattering amplitude in 
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Figure 16: Amplitude Mo (12|12') in Glauber approach. 



this case has the following form 



with 



RA 



(10) 



iVo(lO) = 1 



2 as 



(1 



-il-xio\ 



IF{1) 



(B-1) 



(B-2) 



where / is the transverse momentum of one of two gluons in the Born Approximation diagram. 

The function IF{1) in Eq. ( |B-2| ) denotes so-called impact factor and describes the interaction of two gluon with 
a nucleon which depends on P . In our further calculations we do not need to know the explicit form of IFil). 

The profile function T {b; Ra) determines the number nucleons that can interact with the dipole at a given impact 
parameter b — ^{xi + xq). The profile function T {b; Ra) is defined through the nucleon density in a nucleus p{b, z) 
as 

dzp{b,z); J Sb J dz p{b,z) = A (B-3) 

In the simple model z) — pqQ {Ra ~ b) Q (^\/ R\ — b"^ — z^ , where Ra is the nucleus radius. 

Let us consider a dipole having different sizes at the interaction time t = in the amplitude and conjugate 
amplitude, namely, a dipole with quark (antiquark) coordinates xi {X2) at t = in amplitude and xi {x2i) in the 
conjugate amplitude. A more general case with all the coordinates being different is of no use in our derivation and 
thus will be omitted. The cross section M(12|12') is normalized such that M(12|12) = 2 7V(12). 

As one can see in Fig. |l^ the interaction of the dipole " 12" as well as the dipole " 12"' with a nucleus includes 
an additional elastic rescattering as well as the inelastic interaction with the nucleons of the nucleus at points z^. 

To include both of the processes we first write the formula that sums all possible inelastic interactions. It has 
the form of [El ^ 



M|5" (12|12') = 

(■2R.A(b) 

X 



/ pa™(12|12')dzoe-^('^"'(i^)+-""(i2'))''2^-(^) g / dz,... / dz„_i / 



(B-4) 



r2RA{b) 



C2R.A{b) 



dznP^a^Mm') 



The factor of exp (— i(cr(12) + a {12')) p 2 RA{b)) in Eq. (B-4)describes the fact that there are no inelastic in- 
teractions between points zi, Z2 . . . Zn- Any of the inelastic interactions occurs with the same nucleon in amplitude 
and conjugated amplitude since the contributions of the inelastic processes with different nucleons can be expressed 
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through the nucleon correlations in the nucleus which are suppressed in heavy nuclei. The sum in Eq. (B-4) expo- 
nentiates and gives 



M™(12|12') = (^e+^Um2') T(b;R^) _ i^e-H-"^ 



(12)+£7«'*(12')}T(6;_Ra) 



(B-5) 



a„(12|12') 



1 



-il-x 



X (l - e-^'-^i2') IF{1) 



= - (t7^^(12) + (7^^(12') - (7^^(22')) 



where (22') ^ (t^^(|xi2 - f 12' |) 



Using Eq. ( p3-6| ) and Eq. ( |B-7| ) we can reduce Eq. ( |B-5| ) to 

M™(12|12') = g-i<T«'^(22') T(6;i?.4) _ g-i{<T«-^(12)+a-S'*(12')}T(6;_RA) 



(B-6) 
(B-7) 

(B-8) 



To obtain the final formula we need to add to Eq. (B-S) the contribution of the elastic scattering which has the 
obvious form 

M^' (12|12') = fl - e-i-'-'d^) Tib;R^) \ ^ g-ia--*(i2') T(fc;fl^)\ (B_9) 



Following the lines of Refs. |g, [Iq] we use the explicit expression for the powers of the exponentials in Eq. (B 



and Eq. (B-9), namely, 



\a''^{l2)T{b-RA)^\x\^Ql 
where Qg is the saturation scale used in the McLerran-Venugopalan model ^ ?] defined as 

Ql{h) = 2aspT[b) 



Using this simple notation we can readily find from Eq. ( B-8 ) and Eq. ( B-9 ) 



(B-10) 



(B-11) 



Mo(12|12') = 1 + e-"22'Qa/4 _ ^~xl,Ql/i _ ^~x\^,qIia 
Mo(12|10) = 1 + e-^^oQ^/-^ - e-"'^'^^/^ - e-"ioQ^/4 



Mo(10|10) = 2(1 - e-^i"^'"/"*) 



and other required combinations. 



(B-12) 
(B-13) 
(B-14) 



C. First step of the evolution in the Glauber approach 

In this section we consider the first step of the evolution for the dipole "12" ("12"') having different transverse 
coordinates in the amplitude and the conjugate amplitude, namely, we emit one extra soft gluon "3" as depicted 
in Fig. 0. The evolution in this case is described by the equation Eqs. (11) -(16) and its first iteration gives the 
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required emission of soft gluon "3". We analyze this situation both using Eqs. (4.1) -(4.6) and show its equivalence 
to the direct application of the Glauber theory. 

The soft gluon "3" can be emitted before or after the interaction time t — 0. Firstly, we consider the contribution 
of the diagram A in Fig. |^ where the soft gluon is emitted before r = in both the amplitude and the conjugate 
amplitude. The direct application of the Glauber expressions for the scattering cross section of two dipoles " 13" and 
"32" ("32"') found in Appendix ^ gives in this case 

clastic scattering in the amplitude and conjugated amplitude inelastic scattering in the amplitude and conjugated amplitude 



,)Q'i/'i\ I p-a^^a'Q?/^ _ p-(2^?3+^32+^32')Qa/4 



^1 _ e-(^?3+^32)Q'/4^ ^1 _ e-(^?3+^32') 
= 1 + e^42'Ql/^ - e-(^?3+^32)Q'/4 _ ^-{xIs+xI^,)qV4 



times the Kernel of the dipole splitting (the emission soft gluon "3") given by 



( ^13 


X23\ 














\^13 


^23 / 


\^13 


^2'3 / 



(C-2) 



On the other hand we can use Eq. (12) and Eq. (O) to obtain the same result as follows 



Mo(13|13) + Mo(32|32') + Mo(13|13)A.fo(32|32') - 2Mo(32|32')A^o(13) 

-Afo(13|13) {iVo(32) + A^o(32')} + A^o(13) {A^o(32') + iVo(32)} (C-3) 
^2(1- e-^i3Qa/4^ + + e^"'^"''^'^^ - e-"'32Qs/4 _ g-:i:32''3'/4^ 
2(1- e"^"^'/'') (1 + e""^22'Q'/4 „ e-xl^Q^Ji _ g-:rL'Qe/4^ 



_2 [1 - e""^J3Q;/4 j (^i _ e"=^32Q;/4 _^ ^ _ ^-x'^^,Qi/i^ 
+ (1- e^"^i3Q'/4^ _ e-'^hQl/i ^ ^ _ g-x^,,Q^/4^ 
^ 2 (^1 - e^"^i3Qe/4^ + + e^"^22'Q'/4 _ e-^32Qs/4 _ e-^32'Q'/4^ 

^ 1 + e-^22'Q'/4 _ e-(^?3+^i2)Qs/4 _ ^-{xl;+xl^,)Q^J4 

In the limit of X2 = X2' this equals to the total cross section of the two dipole scattering 2(1 — e~^^i3"''^32)Qs/4) as 
expected from the optical theorem. 

Next, we want to obtain the contribution of diagram B in Fig. |^. The direct application of the Glauber formulae 
gives 

elastic scattering in the amplitude and conjugated amplitude inelastic scattering in the amplitude and conjugated amplitude 



(^1 + e-("'?3+^32)Qa/4^ + e-:r?2'Q'/4^ + e-(''32+^32')Q'/4 _ e-(2^?3+^i2+^32')'3e/4 (C-4) 

= 1 + e"("'32+^32')Q'/4 _ g-(x?3+2;^3)Q2/4 _ ^-^xl^,Ql/4 
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The same result one obtains using Eq. (4.4) in the evolution equation for M(12|12') and the eikonal formulae for 
one dipole scattering 



Mo(43|12') {1 - iVo(32)} + iVo(12')A^o(32) 

= (l + e""=32'Q'/4 _ ^^x^^Q^ji _ g-^?2'Qa/4^ (-^ " ^ + e"^'2'3'/4 



)(> 



-^32Q'/4 



(C-5) 



The contribution from diagram _B* is found in analogous way. 

The remaining diagrams i?, i?*, Z?, E and F change the cross section merely by multiplication of the Kernel 
of the emission of the soft gluon "3", in other words, when gluon "3" is emitted in those diagrams, only the initial 
dipole interacts with the target as 



Mo(12|12') = 1 + e-"'2'«'/4 - e-^l^Qll^ - e-?2'Q^/4 



(C-6) 



In this way we show that evolution equation Eqs. (4.8)-( 4.11 ), correctly describe the evolution to the first order. 
The higher order iterations can be found in the same way. 



D. Glauber expression for single inclusive cross section: AGK cuts 



In this Appendix we show how the Glauber expression for the single inclusive cross section with inelastic part retained 
found in Section]^ (see Eqs. (7.1)-([7^) can be obtained directly applying the Glauber formula. The cross section 
for single inclusive production consists of terms coming from diagram A in Fig. |(see Eq. (O)), diagrams B and 
B* (see Eq. ( |7.4| )) and diagram D (see Eq. (7.4)). We follow the lines of the Appendix ^ and readily write the 
Glauber expression for diagram A in Fig. as 



elastic scattering in the amplitude and conjugated amplitude 

{a-B-*(12)+a-B-*{20)}^ _ ^-i{a-B-4(i2') + <T«'*(2'0)}^ 

inelastic scattering in the amplitude and conjugated amplitude 



(D-1) 



+ e 2 



{<T-^-*(12)+<T^'^(20) + o-^'*(12')+(T^-*(2'0)} -1{ct-^-*{12) + ct-^'*(20)+ct-^-*{120+ct-^"*(2'0)} + ct;„(12|12M+ctj„ (20|2'0) 

— e 2 



_ g-i{<TB'*(12) + £rB'*(20)} 



-i{<T«-*(12')+<T«'*(2'0)} 



- e 2 



_g-i{<T^'*(12)+o--S-*(20) + tT^^(12') + o-S^{2'0)}+o-i„(12|12') + CTi„(20|2'0) 



where we kept 2 7^ 2' (coordinate of gluon in the amplitude and the conjugate amplitude) because we fix the 
momentum of the emitted gluon, and absorbed T{b\RA) in the definition of . In the diagram A in Fig. ^ we 

depicted only splitting, all other kernel have the same Glauber expression. This is not the case as far diagrams 

^12 ^12' 

of the type B, B* and D are concerned, their Glauber expression depend on the dipole splitting, since some inelastic 
contributions are suppressed in large Nq limit. We have already discussed this point deriving the single inclusive 
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cross section Eqs.(3.2)-(3.6) in Section |3| The diagram B, as it is shown in Fig. ^, brings 



clastic scattering in the amplitude and conjugated amplitude 



(D-2) 



inelastic scattering in the amplitude and conjugated amplitude 



+ p-H'^'''^(12)+T-^-*(20)+a-S-4(10)} _ g-i{CT«'*(12)+CT«'*(20)+a«^(10)}+ff.„(12|10) 

^ 1 _ ^-^{<y'''^(12)+a'''^{20)} _ ^-la^^{W) _ g- i {ct^-* (12)+a^'* (20)+a-S-* (10) }+<t,„ (12| 10) 



All other splittings for type B diagrams are found in a similar way by replacing (Tm(12|10) by the proper inelastic 
term, same for B* . The term resulting from diagram D in Fig. is readily calculated as 



clastic scattering in the amplitude and conjugated amplitude 



(l_e-^-""(i"))(l-e-^'^""(i°)) (D-3) 



inelastic scattering in the amplitude and conjugated amplitude 



_^g-i2a«^(10) „ g-i2^«^(10)+*„.(10|10) 
= 1 _ ^-h-^^m _ g-i-^^(lO) _ g-i^ 



The crossed dipole splittings 



and 



in ZJ-type diagrams, has inelastic cross section suppressed in the 
an see that this way we reprod 
in the Glauber form Eqs. (7.1)-(7.4) found from a general expression Eqs. (|3.2|)-(|3. 



X02 ; 

large Nc limit giving just the elastic term. One can see that this way we reproduce the single inclusive cross section 
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